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Abstract. Correlations among the independently measured physical 
properties of globular clusters (GCs) can provide powerful tests for the- 
oretical models and new insights into their dynamics, formation, and 
evolution. We review briefly some of the previous work, and present pre- 
liminary results from a comparative study of GC correlations in the Local 
Group galaxies. The results so far indicate that these diverse GC systems 
follow the same fundamental correlations, suggesting a commonality of 
formative and evolutionary processes which produce them. 



1. Introduction 

Understanding of the physics, formation, and evolution of any type of astro- 
nomical objects or systems, including globular clusters (GCs) as a family, must 
rest on a solid, quantitative empirical foundation. In the order of an increas- 
ing information content, the first step is a definition of characteristic values, 
e.g., a typical mass or luminosity or half-light radius, etc. The next step is a 
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determination of distribution functions for various physical quantities, e.g., the 
luminosity function. Finally, most information can be obtained from non-trivial 
correlations of independently measured or derived quantities. Such correlations 

are products of some physical or evolutionary processes, and as such contain 
valuable clues towards the understanding of the objects in question. 

In the context of globular clusters, previous studies include, e.g., Chernoff & 
Djorgovski (1989), Djorgovski (1991, 1995, 1996), Djorgovski & Mcylan (1994), 
Bcllazzini ct al. (1996), Bellazzini (1998), McLaughlin (2000), and numerous 
papers by Sidney van den Bergh, e.g, van den Bergh (1996). For a comparison 
of GCs with other old stellar systems, sec, e.g., Kormendy (1985) or Djorgovski 
(1993). On-line compilations of Galactic GC data useful for such studies include 
Harris (1996), ^ which is an updated superset of the data ^ presented in the 
First Ivan King Festschrift (eds. Djorgovski Sz Meylan 1993). Both are now well 
overdue for a major revision and updates, including Hipparcos-based distances, 
uniform IR photometry and reddenings derived from the 2MASS data, better 
core parameters from the HST-based surface photometry, etc. 

Additional insights can be gained by comparing correlations of GC prop- 
erties in GC systems of different galaxies. Currently, we are in practice limited 
to the galaxies of the Local Group. Previous studies include, e.g., Fusi Pecci 
et al. (1994), Djorgovski ct al. (1997), Dubath & Grillmair (1997), Dubath ct 
al. (1997), Meylan et al. (2001), Barmby et al. (2002), etc. In this paper we 
present preliminary results from a new study of dynamical correlations for GC 
systems in several Local Group galaxies. 

2. Correlations for Galactic GCs: A Brief Overview 

A generic expectation from the differences in corresponding dynamical time 
scales is that the evolution at the core radius scales would be much faster than 
at the half-light radius scales. The former then reflects mainly the evolution 
towards the core collapse (and any self-similar behavior before and after the core 
collapse), whereas the latter reflects more the initial conditions and a long-term 
evolution. Position of a GC in the Galactic potential modulates the evolutionary 
effect through the effect of dynamical shocks due to the disk and bulge passages, 
in the sense that clusters exposed to more frequent and stronger shocks evolve 
faster. This introduces a secondary dependence on the distance to the Galactic 
center and plane in many of the observed trends and correlations. All this 
was well documented in the references cited above, and at least the qualitative 
agreement between the observations and theory is striking. 

In a simple core collapse picture, as t ^ tc (time of the maximum collapse) , 
the core radius rc — > 0, the concentration parameter c — > oo, the central density 
po ^ oo, and the central surface brightness /q oO; in a self-similar manner; see, 
e.g., Meylan & Heggie (1997) for a review and references. The observed cluster 
properties at the core scale are thus driven towards a (nearly) 1-parameter se- 
quence corresponding to the relative dynamical time away from the core collapse 



^http://physun. physics. mcmaster.ca/Globular. html 
^http://www.£istro.caltech.edu/^george/glob/data.html 
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Figure 1. Examples of core parameter correlations, driven by the 
evolution towards the core collapse. Left : Correlation between two 
observables, the King concentration parameter c, and the core radius 
Tc- Clusters with an unresolved (from the ground) post-core-collapse 
(PCC) morphology have been assigned c = 2.5. Right : An illustra- 
tive correlation between two derived quantities, the central luminosity 
density, po, and the central relaxation time, trc- 



(Fig. 1). This is indeed consistent with the observations: the predicted trend 
in simple single-component, Fokker-Planck core collapse models is pQ ~ r^.^^'^^, 
whereas the observed trends are Iq ~ j,^-i-8±0.2 g^j^^ po^iurn ~ ^^-'2.(>±o.2 ^ This 

implies for the core mass rricore ^ rJ^'^^^''^{M/L), i.e., a nearly constant core 
mass, possibly slightly diminishing due to the evaporation of high-energy stars, 
and/or becoming slightly darker due to the mass segregation of heavy stellar 
remnants. 

No such trends are seen at the half-light radius scale, where the dynamical 
range of the relaxation time scales relative to the Galactic age are is much smaller 
{trh ~ 10^ — 10^'^ yr, whereas for the core relaxation time scales trc ~ 10^ — 10^'^ 
yr), so that the internal spread of GC properties dominates over the dynamical 
evolution effects. 

Likewise, there are only a few noisy trends with the cluster luminosity 
(~ mass), which are apparent only when the data are binned: more luminous 
clusters tend to be more concentrated, ~ ^-0-5±o.25^ ^ L^^^. 

One can also estimate the rough tidal radii rt from the observed surface 
brightness profiles, and the mean cluster densities pt within the rt. One finds 
a mean trend with the present Galactocentric radius Rqc of rt ~ Rq^^"^^'^^, 
as intuitively expected, and also pt ^ Rq^}'^^^ '^ , close to the mean density law 
for the dark halo. Using a simple theory (Innanen et al. 1983) and an assumed 
Galactic rotation curve, it is then possible to estimate the perigalactic radii, 
Rperi (see Djorgovski 1996 for more details). Intriguingly, the ratio Rcc/Rperi 
peaks near the unity, with a long tail, suggesting that most GCs today are on 
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Figure 2. The Fundamental Plane (FP) bivariate correlations for 
GCs: an optimal combination of the central velocity dispersion (cr) 
and core radius (rc) is correlated with the central surface brightness 
in the V band (left panel), and a different combination of a and the 
half-light radius r/j with the mean surface brightness within (right 
panel). Pearson (r) and Spearman (s) correlation coefficients are indi- 
cated in each panel. The residual scatter is completely accounted by 
the measurement errors. 



nearly circular orbits. This could be a survival selection effect, or a reflection of 
the initial conditions, or a combination of the two. 

Possibly the most interesting monovariate correlations are those including 
core velocity dispersions, cr: cr ~ ^o.6±o.i5 ^^j. equivalently, L ~ cr^/^), a ~ 
j^o.5±o.i^ and a ~ j^o.45±o.05^ where 1^ is the mean surface brightness within 
the half-light radius (see Figs. 3-4). The origin of these correlations is still not 
understood, but several possibilities exist; see, e.g., Djorgovski & Meylan (1994) 
for a discussion. These correlations may be reflecting the formative mechanisms 
of GCs, and they should be reproduced by any successful model of GC formation 
and long-term evolution. They are also very different from the corresponding 
correlations for elliptical and dwarf (dE, DSph) galaxies (Djorgovski 1993). 

It is notable that GC metallicities do not seem to correlate with anything, 
in contrast to other old stellar systems. This strongly suggests that GCs were 
not self-enriched. 

Another interesting question is how many independent physical parameters 
control the observable properties of GCs? Following the pioneering study by 
Brosche Sz Lentes (1984), more modern data sets suggest that there are at most 
6 significant parameters among 9 or 10 independently measured observablcs 
(Djorgovski 1991, Djorgovski &: Meylan 1994). If one considers just the observed 
photometric, structural, and dynamical parameters at both core and half-light 
scales, the statistical dimensionality of the data is 3, or 4 if the {M/L) ratios 
are included. This is exactly as expected for a manifold of King (1966) models. 
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Figure 3. Luminosity - velocity dispersion correlations in the V (top) 
and K (bottom) bands for GCs in 4 Local Group galaxies, encoded 
with different symbols as indicated in the figure. The clusters in N205 
appear to be slightly underluminous for their velocity dispersions, but 
this may be due to a systematic error. Otherwise, all of these GC 
systems appear to follow the same dynamical correlations. 
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Figure 4. Correlations between the velocity dispersion a and the 
mean surface brightness within the half-light radius in the V band 
(HhVi Isft panel) and in the K band {^hK-, right panel) for the Galactic 
GCs (open circles) and M31 GCs (solid squares). The M31 extends 
to higher luminosities and masses, and is selection-limited at the low 
end, but in the overlap region the two GC systems follow the same 
correlations. 
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Constraining the input parameter set to cither core or half-light ones, brings 
the statistical dimension of the data set to 2, i.e., the triplets of observables are 
connected by bivariate correlations (Djorgovski 1995; see Fig. 2). This is the 
equivalent of the Fundamental Plane (FP) of elliptical galaxies. For the core 
parameters, the scaling relation is: Tc ~ ^2.o±o.2j^-i.i±o.i^ j£ assume a 
structural homology (probably a good assumption for the GC cores), the virial 
theorem implies: r ~ a'^ I^^ (M / L)~^ . Thus, the FP of GC cores implies that 
their (M/L) ratios are (nearly?) constant; this is the sharpest observational 
constraint on the constancy of GC (M/L) ratios to date. In contrast, for the 
half-light parameters, the scaling relation is: r/j ~ ^i.45±o.2j^-o.85±o.i_ 'pj^jg jg 
remarkably similar to the FP of E- galaxies, and is almost certainly a consequence 
of the non-homology of GC structures. 

An alternative look at the FP of GCs was provided by McLaughlin (2000) . 
While recognizing the equivalence of his FP with the one described above, 
his preferred scaling relation is the expression of the binding energy, Ef) ~ 
L^ °^i?^'*, where the 2nd term effectively corrects for the known dependences 
of the GC parameters on their position in the Galaxy. However, wc note that 
the binding energy can be written as: Ej^ = L!^ [M / L^r f {c) , where /(c) is 
a slow function of the cluster concentration. If the {M/L) ratios are indeed 
constant, and knowing that Vh does not correlate much with anything, it is then 
not surprising that the observed i?GC-corrected scaling is E}, ~ L^ °^. 
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Figure 5. Dependence of the derived (M/L) ratios (in Solar V band 
units) on the metallicity, for the Galactic GCs (open circles) and M31 
GCs (solid squares), in the V (left) and K (right) bands. The larger 
observed scatter for the Galactic GCs may be caused in part by the 
distance errors, whereas all M31 GCs are effectively at the same dis- 
tance. While no trend is seen in the V band (as implied by the FP; see 
above), there is a trend in the K band in the sense of more metal-rich 
clusters having more luminous stellar populations at a given mass. 



3. A Comparison of Dynamical Correlations for GCs in Five Local 
Group Galaxies: Some Preliminary Results 

Correlations between the velocity dispersion and other parameters (L, Iq, I^, and 
the FP) may probe directly the physics and formative processes of GCs, their 
homogeneity (or lack thereof) in different galaxies, and can be also used as new 
distance indicator relations for their parent galaxies, providing an independent 
check of other distance scales. Fusi Pecci et al. (1994) and Barmby et al. (2002) 
established that the structural and photometric parameters of GCs in M31 as 
observed with the HST occupy the same portion of the parameter space as the 
Milky Way GCs. Following on the initial studies by Dubath & Grillmair (1997) 
and Djorgovski et al. (1997), we set to explore in more detail the dynamical 
correlations for GCs in 5 Local Group galaxies: the Milky Way, M31, M33, 
N185, and N205. 

The extra-Galactic sample consists of HST imaging of GCs in M31, M33, 

N185, and N205 available as of the early 2002, for which we have done sur- 
face photometry (Federici, Parmeggiani, et al., in prep.). Velocity dispersions 
have been measured using Echelle spectra obtained at the Keck-I telescope with 
the HIRES instrument, as described, e.g., in Djorgovski et al. (1997); these 
measurements will be presented in detail elsewhere (Cote et al., in prep.). 

Figures 3-5 illustrate some of our preliminary results. Our first conclusion 
is that GCs in these different galaxies follow essentially the same dynamical 
correlations. This suggests a common set of physical mechanisms affecting their 
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formation and evolution, despite a broad range of their host galaxy properties. 
The observed dependence of the (M/L) ratios on the metallicity, especially in 
the K-band (Fig. 5), represents a useful observational constraint on the models 
of old stellar populations. In a future paper (Djorgovski et al., in prep.) we will 
present a more complete and detailed analysis, including the FP correlations for 
these GC systems. 

Even with the resolution of the HST, such studies cannot be pushed much 
beyond the Local Group. For example, Harris et al. (2002) present an excellent 
study of GCs in Gen A = N5128. 
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